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1 Introduction. 

Very soon after publication of the famous BCS work explaining the puzzle 
of the superconductivity, N.N. Bogolyubov and coworkers JIJ have proposed 
their version of the theory. One of new results they have obtained was the 
discovery of a collective mode - an oscillation of the Cooper pair density with 
the energy smaller than 2 A. P.W. Anderson has indicated that this collective 
mode can not be found experimentally since the Coulomb forces neglected 
in the above mentioned work shift its energy to the plasma frequency, i.e., 
to the high-ultraviolet range in which the superconductivity is inessential. 
The interest for plasmons in superconductors revived a little in the search 
for the mechanism of the High-T c superconductivity. However, a real surge 
of interest to this problem occured after experimental observations of the 
plasma edge in the reflectivity of High-T c superconductors La 2 ^ x Sr x CuO^ 
and YBa2Cu 3 8 - y @, 0. It is worthwile to mention that a theoretical 
prediction have preceded the experiment |4j] . For external reasons, the paper 
M has been published only a long time after its completion. 

Here we present a brief review of the experiments and theoretical devel- 
opments in the field. The theoretical works will be presented in more detail, 
given the author's specialization. I hope that this imbalance will be compen- 
sated by the papers by O.K.C. Tsui, N.P. Ong and J.B. Peterson, and by 
P.M. Mueller. 

The content of the review is as follows. In the second section we review the 
most general relationships for plasma frequency (PF) and dielectric function 
of the electron plasma in normal metals and their modification for layered 
normal metals. In the third section we review the experimental observations 
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of the plasma edge in High-T c superconductors (HTSC). In the fourth section 
a simple two-fluid model of the plasma oscillations is presented. The fifth 
section is devoted to the BCS theory of the same phenomenon. Peculiarities 
of Josephson-linked layered superconductors are considered in Section 6. The 
experimental observations of the resonance plasma attenuation are discussed 
in Section 7. Mechanisms of dissipation in the Josephson-linked supercon- 
ductors are discussed in Section 8. A novel effect of selectional transparency 
is discussed in Section 9. The conclusions are collected in Section 10. 
We employ abbreviations SC for superconductors and OP for the order pa- 
rameter. 



2 Plasma resonance in normal metals. 

The plasma oscillation is one of the simplest and most fundamental manifes- 
tations of the Coulomb interaction. Due to the long-range character of the 
Coulomb forces, its spectrum starts with a rather high frequency given by 
the well-known formula: 

nj = — (i) 

where n and m are the electron density and mass respectively, is the 
dielectric constant contributed by atomic and ion reminder shells. For usual 
values in solids n ~ 10 22 -j- 10 23 cm -3 , m = 10~ 27 g and ~ 1, the PF 
Q p i ~ 10 15 -J- 10 16 s _1 lies in the high-ultraviolet range. The variation of n 
and m in metals are not large. Therefore the characteristic value of the PF 
is always rather high. This is the reason why plasma oscillations do not play 
an important role, either in the thermodynamics or in transport phenomena 
in conventional metals. 

The situation changes in layered metals, if the electron propagation per- 
pendicular to the layers (the c-axis) is mediated by tunneling processes. Ap- 
plying the tight-binding approximation for the c-direction one finds a follow- 
ing dispersion law for electrons: 

<P,Pz) = e||(p)+*(l-cos^) (2) 

where p,p z are in-plane and perpendicular-to-plane components of the mo- 
mentum, t is the overlap integral and eii(p) is the in-plane dispersion. As- 
suming it to be eii(p) = p 2 /2m for simplicity, we find that for t < ej? the 
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Fermi-surface is a slightly corrugated circular cylinder. Solving the kinetic 
equation with such a Fermi-surface, we find that Eq. ([I]) holds for the electric 
field in-plane, but for the perpendicular-to-plane electric field the PF is much 
smaller j5|: 

^ = (3) 



where M = h 2 /td 2 m and d is the lattice constant in the c-direction. 



Note that the ratio Q c /Q a b oc m/M instead of ym/M as it could be naively 
expected and what is correct for an elongated ellipsoid. The reason is the 
change of topology: the cylinder is an open surface. This simple observation 
allows one to understand a drastic reduction of the PF which has been ob- 
served experimentally 0, p): Q p i for LaSrCuO and for YBCO is in the far 
infrared region. This agrees well with a comparatively low carriers densities 
found in the Hall-effect measurements and high ratio M/m ~ 100 -j- 2,500 
found in experiments with the tork produced by the vortex lattices 0. The 
anisotropy parameter in Bi 2:2:1:2 is even larger ||, M/m ~ 40, 000. There- 
fore we predicted that plasma frequency in this compound lies in the mi- 
crovawe region ||. This theoretical prediction has been confirmed in the 
experiment ||, though the experimenters were not aware of the theoretical 
prediction. 

Unless the collision time for quasiparticles r is small enough, there is no 
substantial difference between the plasma oscillations in normal metals and 
superconductors. In both cases they are the charge density oscillations with 
the restoring force determined by the electrostatics. The situation is very 
different for small collision time Q p it <C 1. Then normal carriers almost do 
not participate in the plasma oscillations. We return to the analysis of this 
problem later. 

The PF may be alternatively determined as a frequency at which the 
dielectric function e(u) is zero. For the collisionless plasma, the dielectric 
function is: 

e(.) = , x (l - §) (4) 
The reflectivity for an EM wave is determined by the Fresnel formula: 
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e(u) + 1 



(5) 
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According to Eqs. (||, |^) R = 1 at uj < Q p i and R < 1 at uj > Q p i. This 
effect is called the plasma edge. Its sharpness depends on the value e^. For 
> 1, the graph of e(u) has a characteristic beak (Fig 1). Introducing a 
phenomenological collision time r and the electron attenuation V = 1/r, one 
finds instead of (H): 

The plasma edge is smeared out at large T 3> Q pi or small r. Theoretical 
curves of R(u) are shown on Fig.l. 

3 Experimental observation of the plasma edge 
in superconductors. 

The first observation of the plasma edge has been reported by Koch et al. 
[0] in 1990. They used a single crystal of YBCO specially grown in such a 
way that the layers were oriented perpendicularly to the surface of a sample 
and the c-axis was directed parallel to the surface. The electric vector in the 
incident EM wave was directed along the the c-axis. This geometry enabled 
the authors to observe the plasma edge. They have found a sharp decrease 
of the reflectivity at a frequency uj ~ 12meV in the superconducting state 
(T < 90K) and very smooth decrease of reflectivity in the normal state (Fig. 
2). Unfortunately, the plasma edge in YBCO occured to be close to other 
features of the spectrum associated with optic phonons, thus interfering with 
a thorough study of plasma oscillations. 

Tamasaku et al. || have used very good single crystals of La2- x Sr x CuO^ 
in the same geometry which enabled the authors to observe the plasma edge 
reliably (Fig. 3). They also did not observe any edge effect in the nor- 
mal state and a rather sharp edge already at a temperature just slightly 
below the transition temperature. A systematic decrease of the PF with the 
temperature has been found. Their measurements are compatible with the 
assumption that oc A 2 . 

Two questions naturally arise in connection with these measurements. 
1. Why there is no plasma edge in the normal state while it occurs immedi- 
ately in the superconducting state? 
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2. Why PF depends strongly on the temperature? 

At a qualitative level these questions can be answered without detailed cal- 
culations 0, |J. The absence of the plasma edge in the normal state can 
be explained by strong collisions VL p iT C 1. In the superconducting state a 
part of electrons is confined in Cooper pairs (superconducting carriers) which 
are not scattered. It explains why the plasma oscillations occur in the su- 
perconducting state. The same idea explains the temperature dependence 
of the PF: not the full, but only the superffuid density n s enters Eqs. ([[], 
H). The superffuid density in an anisotropic SC is a second-rank tensor. The 
superffuid density varies from the total electron density at T = in a clean 
SC till zero at T = T c . 



4 Two-fluid theory of the plasma resonance. 

The two-fluid model of SC by Gorter and Kasimir (see for example is a 
simplest phenomenological way to understand the plasma oscillations in SC. 
It has been used in fact by T. Mishonov j|, and also in a great extent by D. 



van der Marel et al. [II] and later by Tachiki et al. [12]. In the framework 



of this model, the dielectric function is written in the following form: 

eM-eoofl-gf- , ";- n V (7) 
\ or to(u + iF) J 

The terms with and Vt^ represent contributions of the superffuid and 
normal parts of the electron liquid, respectively. As it was emphasized in 
Section 3, the superffuid part of the liquid is not subject to any dissipation 
whereas the normal part contains the conventional dissipative coefficient V. 
If fl p iT ^> 1, then Q pl — Q% + fl^, i.e. all carriers participate in the plasma 
oscillation. In the opposite case Q p it <C 1, one finds fl p i = Q s , in agreement 
with what was expected. Note that in the latter case the imaginary part of 
e is small: 

Jme = !^I. {uj *n s ). (8) 

Thus, the stronger are collisions the weaker is the dissipation of the plasma 
oscillations. This paradoxical result is explained by different roles played 
by normal carriers at different Q p it. If this parameter is large, the normal 
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carriers create the coherent plasma wave on the same footing as the super- 
conducting electrons. If Q p it is small, the coherent motion is created by 
superconducting electrons only whereas the normal electrons, being partly 
involved into this motion, dissipate the energy. The sooner they are put in 
equilibrium by collisions the more adiabatically they follow an instant value 
of the OP, and the smaller are losses. 

Though the two-fluid model describes qualitatively the plasma oscilla- 
tions it has several important failures. First, the very notions of the su- 
perconducting and normal densities n s and n n can be defined reasonably 
only for a rather clean SC. These definitions presume the conservation of the 
momentum which is violated by impurities, phonons etc. Second, the phe- 
nomenological definition of the collision time does not enable one to calculate 
its temperature and frequency dependence. Third, the two-fluid model does 
not take into account the symmetry of the superconducting OP making no 
difference between s- and <i-pairing. 



5 BCS theory of plasma oscillations. 

The BCS theory of plasma oscillations has been proposed in || to overcome 
the failures of the two-fluid model. It is a microscopic theory based on the 
BCS Hamiltonian with the dispersion law given by Eq. (0). The OP in 
the unperturbed SC is assumed to be uniform in the real space and very 
anisotropic in the relative momentum space. Only the elastic scattering by 
impurities has been considered. The scatterers were assumed to change the 
electron momentum in plane, but to be unable to transfer an electron from 
one plane to an ajacent plane. It means that the scattering amplitude does 
not depend on the z-component p z of the momentum. One more simplifica- 
tion is provided by an assumption of the full in-plane isotropy. The problem 
is to calculate the Pippard kernel Q(u,q). The Pippard kernel relates the 
Fourier components of the current density j and the vector potential A: 

j(w,q) = -Q(w, q) A(w, q) (9) 

The relationship is written for a special gauge VA = 0. Once the Pippard 
kernel is known the dielectric function can be calculated as: 

Aire 

e(w,q) = -Q(cj,q) (10) 
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In the real space the relationship (|9[) reads: 

j a {w,r) = - J Q a)(3 (uj,r,v')A/ 3 (uj,r')d 3 x 



(11) 



The relationship (JTTJ) is valid at a fixed configuration of impurities which 
violates the translational invariance. According to the Kubo formula, the 
response function (Pippard kernel) is expressed in terms of the retarded 
current- current Green function which in turn can be obtained as an ana- 
lytical continuation of the Matsubara Green function: 



Qa 



M 



,UJ r , 



iT 



l/T 



dre^(j a (r,v)j p (0,v')) 



(12) 



with u n = 2irnT and integer n. The retarded response is the analytical 
continuation of this function from discrete points uj = iu n in the upper half- 
plane of the complex plane u. The next averaging over the random location 
of impurities should be performed in the end. 

Abrikosov and Gor'kov [DJ have developed a special diagrammatic tech- 
nic for this averaging. In a simplest situation of weak scatterers, the problem 
is reduced to the solution of an analog of the Boltzmann kinetic equation for 
SC [14 1, [HJ. An essential complication for the superconducting case is the 



existence of several distribution functions instead of one in the normal case. 
Not only the occupation numbers of quasiparticles, but also the condensate 
wave function varies locally. The solution of the kinetic equation is gener- 
ally nontrivial even in the normal case. Fortunately, our special geometry 
together with the fact that the penetration depth is much larger than any 
other geometrical scale in the problem makes this problem exactly solvable. 
More accurately, it can be solved in the leading order over small tunneling 
amplitude t if one employs the independence of the impurity scattering am- 
plitude on p z . For the most interesting zz-component of the Pippard kernel 
the result is: 



Q,z(n,T) 



AncTi 



(13) 



+ tanh 



dr\ 

v 



i + 



r] + r] + A 2 \ tanh(r/+/2T) - tanh(r//2T) 



e R (rj + )e A (rj) , 

r] + ri+A 2 

<^r{v+)^r{v) 



e A (rj) +i/r 

r}-r)+A 2 



+ 



eA{rj-)e A (rt) 



2T \ e R (7] + ) + e R (r}) +i/r e A (r}-) + e A (7]) -i/r 



(14) 
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where rj± = 77 ± u, e^A = J (tj ± iS) 2 — A 2 , the branch of the square root is 
chosen to be positive at 77 > A. The reader is referred to the original work 
U for details. 

One of the most important results in the BCS theory of the plasma oscil- 
lations is the vanishing of the dissipation in a pure SC at any T < T c . Two 
conventional dissipation processes are the Cherenkov absorption of light by 
quasiparticles and the Cooper pair breaking, i.e., the creation of a pair of 
quasiparticles with the absorption of a photon. The first process is kine- 
matically forbidden since the light velocity c/eoo is much larger than the 
Fermi velocity. The second process, specific for SC, is kinematically allowed 
at uj > 2 A. However, it is forbidden dynamically. The amplitude of the 
quasiparticle-photon interaction is proportional to the so-called Bogolyubov 
coherence factor |L0| (u p+(l v p — v p+cl u p ) 2 where p is the electron momentum 
and q is the photon pomentum. For the extreme of the normal skin-effect 
q — > this factor is zero. Both processes are allowed in the impure SC in 
which the distribution over momentum at a fixed energy looses its delta- 
like character and becomes a Lorenzian one. This argument explains why 
the plasma edge is observable in the HTSC which have been experimentally 
proven to have lines of nodes of the OP and zero energy gap in the spectrum 
0,0,0,0,0. 

Eqs. ([TU], |TJ]) show that what is called fl s and Q n in the two-fluid model 
in the microscopic BCS theory are complicated functions of all variables. 
Even at very low temperature the analog of the superfluid density n s depends 
rather sharply on r varying from the value of total electron density at r A 3> 1 
to a value less by a factor VtA at rA <C 1. 

If Q p i <C A, the PF is expressed in terms of the penetration depth: Q p i = 
c/X c . This is just the case for Bi 2:2:1:2. However, in YBCO and LaSrCuO, 
fl p i is comparable with A. In this case, the PF can be found as a solution of 
a following equation: 

1lee(tt pl ) = Q (15) 

For this situation the exact result (|14j) is especially important. 

In a clean SC, ImQ = 0, as it has been shown earlier. Therefore, K — 1 at 
t = 00, independently on temperature and frequency. This is a generalization 
of the relationship Q pl = Q 2 + f2 2 at T = in the two-fluid model. The 
transition between the static value K oc n s (T) and the pure limit K — 1 
occurs in a region of frequency u ~ 1/r, very narrow in clean SC. 
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Comparing BCS theory and the two-fluid theory of EM response in lay- 
ered SC we find that quantitatively the two-fluid model is rather unsatisfac- 
tory. It is especially poor in evaluation of the dissipation for which its error 
(a deviation from the BCS theory) may be 1000% and more. The reason is 
that there is only one scale of frequency 1/r for dissipation in the two-fluid 
model whereas the BCS theory displays several equally important scales: 
A,T,l/r,A\(Ar 2 r\ 

Though the BCS theory is a step forward in our understanding of the 
phenomenon its premises are oversimplified. The most important of these 
simplifications are as follows: 

i) The BCS theory is a theory of weakly interacting electrons. It may not work 
well for HTSC. Unfortunately the existing theories of strongly correlated 
fermions do not reach a level wich allows for a reliable calculation of the EM 
response. 

ii) The scattering in the BCS model has been assumed to be elastic scattering 
by weak impurity scatterers. It probably dominates in the LSCO which 
becomes a SC at a sufficiently high doping only. However, measurements of 
the Hall angle in YBCO at low temperatures have shown a high purity 
of this compound. In combination with the plasma edge observation this 
leads to the conclusion that the scattering rate is large near the transition 
temperature (Q p it <C 1) and drops fast with decreasing the temperature. 

The problem of the scattering mechanism stays in a close connection 
to the problem of the strong electron interaction. The famous linear- in- 
temperature resistivity p2] can hardly be explained in the frameworks of the 



Fermi-liquid theory. For our purpose we need the c-axis resistivity which be- 



havior resembles semiconductors more than the ||23|| . We can argue that the 
variable-distance hopping processes determining the static c-axis resistivity 
are irrelevant at plasma frequency since they require a long time. Neverthe- 
less, the direct theoretical approach is much desirable. 

iii) Our version of the BCS theory does not take into account a strong spatial 
modulation of the superconducting OP. However, this is rather important for 
explanation of a strong influence produced by the magnetic field on the PF. 
This question is discussed in the next section. 

iv) We have already mentioned that the experiments indicate a nontrivial 
character of the pairing in HTSC with the nodes of the OP. The BCS theory 
should be reformulated for this situatuion. 
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6 Plasma oscillations in Josephson-linked SC. 



In real HTSC the OP is deeply modulated in the c- direction with the lattice 
periodicity. The coherence length in the c-direction, £ c , for Bi, Tl and Hg 
based SC is much less than the lattice constant d in the same direction 
practically at any T < T c . For YBCO, £ c becomes smaller than d at T « 
75K while T c = 90K. The OP in deeply modulated superconductors is 
concentrated in thin layers near conducting planes. The coupling between 
layers is a weak Josephson link. The modulus \ip n \ of the OP in n-th plane is 
a rigid variable wich fluctuates only weakly except very close to the transition 
temperature, whereas the the phases (p n are Goldstone variables, which are 
strongly fluctuating and are easily influenced by external fields. In particular, 
the phase variables are drastically influenced by the magnetic field. This is a 
general mechanism for a strong dependence of the PF on the magnetic field. 

Specifically, the magnetic field penetrates into a layered SC as an Abrikosov 
vortex lattice. We start with a geometry in which the magnetic induction is 
parallel to the layers. Then the vortices form a 2-dimensional lattice with 
central lines situated just in the middle between two ajacent planes. These 
vortices are similar to the Josephson vortices occuring in a Josephson junction 
[E3J]. Each of them creates an additional phase difference between ajacent 
layers 27r on a characteristic scale Aj = 7<i, the so called Josephson screening 



length (JSL), where 7 = yM/m = £ a b/£ c is the anisotropy coefficient. The 
latter is about 30 4- 50 for YBCO and not less than 200 for Bi 2:2:1:2, as 
it was established by several independent experiments. Thus the JSL for 
BiSCCO reaches a macroscopic value Aj ~ 3000A. If the magnetic field is 
so strong that the vortices strongly overlap, i.e., the distance between their 
centers is much less than Aj, the phase in the plane grows almost linearly 



with the distance ||25|| . In the opposite limiting case of weakly overlapping 
vortices, the phase difference Ay? njn+ i between n-th and n + l-th layers grows 
rapidly by 2ir on the scale of JSL near each vortex and remains a constant 
until the next vortex line. The boundary between these two regimes is the 
value of magnetic field Bj = o /7<i 2 . 

To understand how the magnetic field influences the PF we first estab- 
lish an important relationship between the PF and the maximal Josephson 
current j max . The analog of the Pippard kernel Q zz or n s in the Josephson 
linked SC is a value proportional to the maximum Josephson current j ma x- 
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Indeed from the Josephson relationship: 

jz = jmax sm(A(p n>n+1 - ^4^) ( 16 ) 

<Po 

in a special gauge Ay9 nin+ i = and at small A z we discover the linear de- 
pendence: 

jz = ~QzzA z ] Qzz = 7 • (1^) 

(po 

By analogy with Eq. ( |10]) we find: 

Q 2 = 87TC Jmax d 
Coo 00 

The Josephson current can be found by differentiating the corresponding free 
energy : 



2cd 



1 (a 2A ^ 

1 - COS £Hp n>n +l 7 



(19) 



A regular lattice of vortices stretched along the y-axis causes a regular vari- 
ation of A(p nn+ i with x. We consider first B ^> Bj. Then 



A^, n+1 « 2^ ^ + § sin A x j (20 ) 

After averaging of Eq. (|19"D over x with the A<yj n ri+1 given by Eq. (^0]) we 
find 

Je// = jmax-^ (21) 

where j e // is the effective maximal current after the averaging. In the oppo- 



site limiting case B <C Bj the in-plane distance between vortices is J^^o/B 



the out-of-plane distance between vortices is y0 o /7-B and the range of es- 
sential variation of phase near each vortex line is d\j. Therefore 

Jeff=Jma X (l-0(B/Bj)). (22) 

One should expect a sharp dependence of the PF on the direction of the 
magnetic field. Such a dependence has been predicted in pB[ . In particular 
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the authors of the work |26[ have predicted a reentrant behavior of the PF vs 
the angle similar to the reentrant behavior of the ac magnetic susceptibility 
and resistivity discovered earlier P7|, [EH]. 



Next we consider the magnetic field directed perpendicularly to the layers. 
Each vortex in a layered SC consists of separate plane vortices - pancakes. 
If there is no irregularitiy in the location of pancakes, the phase difference 
^■fn,n+i created by them is zero and the magnetic field does not affect the PF 
until B <C H C 2. However, any disorder in the pancake arrangement creates a 
random phase difference and reduces efffectively j max and Q p i. We consider 
the case of the strong magnetic field when many pancakes enter a circle of 
the radius Aj f29[| : 

= & < 23 > 



The condition ( E3|) is much less restrictive than the analogous condition for 
the parallel field: B J± ta 200 Gs for Bi 2:2:1:2. If the condition (|2S|) is 
satisfied, the phase A<p„ in+ i is effectively a sum of a large number N = 
B/Bj± of random terms. Therefore it obeys the Gaussian statistics. The 
effective maximal current is proportional to the Debye- Waller factor: 

jeff = jmaxexp (- - ( ( A^ n , n+ i) 2 )^j . (24) 

The value A<p ni „ + i in our case is: 

A<p n>n+ i(r) = J2 [ $ ( r ~ r n,u) - $(r - r n+ i >v )} 

V 

= J (p n (r)-p„ +1 (r))$(r)rf 2 a; (25) 

where $(r) = arctan(y/x) and r n}U is the vector coordinates of the center of 
a z/-th vortex in the n-th plane; p n (r) = J2u ^( r — r n,u) is the vortex density 
in the n-th plane. After the averaging we find: 

(A^, n+1 ) = 2uln(A J /a); // = -J r 2 /C(r)rf 2 x (26) 

where K.(r) is a correlator: 

/C(r) = (p„(r)p n (0)) - (p„(r)p n+1 (0)> (27) 
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If there is no correlation between ajacent layers, the correlator /C(r) turns into 
the correlator of a two-dimensional pancake liquid. Thus, the measurement 
of the field dependence of the PF provides an important information on the 
statistical properties of the pancake liquid. Assuming the disorder within 
a plane to be so strong that the only characteristic scale for /C(r) is the 
average distance between pancakes a, we find that the dimensionless value \i 
is a numerical constant. As a consequence we find a power law decrease of 
the PF with the magnetic field: 

flj = Q 2 pl0 (B J± /Br (28) 

A reason for a strong in-plane disorder may be either thermal fluctuations 
in the pancake liquid or a rather strong pinning in a glassy phase. The first 
case corresponds to a temperature above the irreversibility line, the second 
one - below this line. In a simplest irreversible situation a large number of 
the strong pinning centers provide a random location of the pancakes. A 
metastable state depends on the way of its preparation. If the magnetic field 
is frozen, the vortices find closest pinning centers and do not deviate much, 
at least locally, from the regular Abrikosov lattice. If the field switches on 
after the cooling, the vortices enter consequently at the field is increased. 
The resulting vortex arrangement may be highly chaotic. 

Consider a situation when there is no correlation between positions of 
pancakes in ajacent planes, but the in-plane crystal order persists over a 
distance R much larger than the in-plane lattice constant a. We also assume 
that R <C A j. Then the exponent in Eq. ( p8|) (see also Eq.(^6|)) becomes 
very sensitive to the model of disorder and may be either much larger or 
much smaller than one, as well as of the order of one. As a consequence, 
the reduction of the PF may be much stronger than that caused by a strong 
in-plane disorder. This striking result is explained by the accumulation of a 
random interplane phase shift on a large area cx R 2 . 

Thus the theory predicts a power-law decrease of the PF at B > B J± 
and a strong in-plane disorder with the exponent of the order of unity. The 
same result may be correct in the absence of the interplane correlation, but 
I can not exclude a possibility of the strong reduction of the PF in this case. 
Finally, one should expect muvh weaker magnetic field effect, if the interplane 
correlation is strong. 

The description of the Josephson-linked SC accepted in this section is 
a kind of a long-wave-length phenomenology basically the same as in the 
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Lawrence-Doniach model |J0|. A microscopic approach has been proposed 
by Artemenko and Kobel'kov [31 . 



7 Experimental observation of the resonance 
attenuation at the PF. 

The first observation of the resonance attenuation has been reported by O.K. 
Tsui et al. ||. A sample of Bi2Sr2CaCu20g+s (Bi 2:2:1:2) has been placed 
into a waveguide. The sample was grown in the conventional way with layers 
parallel to its surface. The resonance has been observed in a range of frequen- 
cies between 30 and 50GHz and magnetic fields from 1 to 7T. The absorption 
was recorded by a resonance bolometric technique. During each experiment, 
the frequency was fixed and the magnetic field varied. Rather sharp lines 
of attenuation have been observed in the temperature interval between 2.5 
and 20K. The effect has been initially ascribed to the cyclotron resonance. 
However, an additional study has shown that it occurs only if the electric 
field has a non-zero c-component, and tha the resonance frequency decreases 



with the magnetic field [ 32| , [ 33j . This fact can be considered as a convincing 
evidence in favor of the plasma-resonance origin of the phenomenon. 

Accurate measurements of the resonance |32|] may be summarized by an 



empirical formula valid at T below the irreversibility line: 

n 2 pl (B,T) = AB-»exp(T/T ) (29) 

where A is a constant, /i ~ 0.8 4- 1.0 and T ~ 12.5K. In Section 6 we 
have derived the power low for at a condition of strong in-plane disor- 
der. Surprisingly the PF grows with the temperature whereas according to 
a naive idea it should decrease as the superfluid density. This growth can be 
qualitatively explained in terms of partial thermal depinning of the vortices. 
Depinned vortices tend to restore the positional order destroyed by the pin- 
ning. Quantitatively we were able to prove that the linear in T correction to 
the PF is positive [29] . 



In Bi 2:2:1:2 the contribution of quasiparticles to the PF is probably 
always small. The reason is that they are strongly damped at T close to T c 
whereas the normal density is very small at T <C T c . Besides of that Q p i -C A 
in Bi 2:2:1:2. 
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8 Dissipation mechanisms in Josephson-linked 
SC. 



The resonance attenuation at lu = Q p i requires an explanation. A most 
obvious effect is the amplification of the electric field inside the SC. In the 
experimental geometry the sample surface was parallel to the layers, the elec- 
tric field in the wave-guide was perpendicular to the surface and the layers. 
From the continuity of the component of the electric induction D z normal 
to the surface, it follows that E int = E ext /e(uj), where E int is the field in- 
side the sample and E ext is the field in the wave-guide. Since e(u) is almost 
zero at u = Q p i, the electric field is strongly amplified inside the sample. 
It creates a necessary premise for the resonance attenuation, but this is not 
sufficient. The problem is that the electric field in the z-direction creates 
the ac Josephson current in the same direction, which is not subject to dissi- 
pation. A mechanism transforming the electric field perpendicular to layers 
into the parallel one and the dissipation mechanism for the parallel compo- 
nent of the electric field or current should be found. Though this question is 
not fully solved we present here preliminary ideas due to L. Bulaevsky and 
S. Pokrovsky without detailed calculations. 

Strongly pinned randomly located pancakes create large fluctuations of 
the dielectric constant e zz . The amplitude of this fluctuations can be es- 
timated as l/(cos Ay9 nn+1 ), the inverse Debye- Waller factor calculated in 
Section 6. The fluctuating dielectric constant leads to a strong and random 



refraction which represents just the necessary transformation mechanism [|3£ 



As for dissipation mechanism two ideas are discussed. The first [39] is the 



dissipation by the normal component which exists in SC with the nodes of 



the OP in magnetic field even at zero temperature [34]: the normal carriers 



appear near vortices along special directions. The density of normal carriers 



is estimated as nJ B/H C 2- This mechanism leads to a linewidth Au/lu ~ 0.1 
in agreement with the experiment. However, it is not clear whether the 
growth of the relative width of the resonance line with decreasing temperature 
predicted by this mechanism is compatible with the experimental data M, 



The second dissipation mechanism is via the excitations of the phase os- 
cillations in superconductors with disordered vortex arrangement. As it was 
discussed above the disordered vortices generate strong static phase fluc- 
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tuations. Due to nonlinearity of the Josephson energy, these fluctuations 
may create inhomogeneous localized phase oscillations. They are the same 
plasma oscillation in an inhomogeneous situation. One can expect a contin- 
uos spectrum of these oscillations starting with zero energy. The excitation 
of these oscillations represents a kind of the inhomogeneous line broadening. 
Estimates made by L.N Bulaevsky et al., give the relative line-width 



~ 0.1. They also predict that the relative line- width does not depend on the 
temperature. 

In the same preprint, Bulaevsky et al. have also considered the well- 
known Bardeen-Stephen mechanism |35| applied to the vortices oscillating 
near their equilibrium pinned positions. It should be taken into account that 
this dissipation is reduced by the quantization of quasiparticles in the vortex 
core. The quantization is insubstantial for low-temperature SC with the core 
size £ ~ 1,000A However, in HTSC ~ 30A and the characteristic level 
spacing is about 100K. Since the experiments have been performed starting 
from 2K the excitation of quasiparticles bound to vortex cores seems to be 
irrelevant. 



9 Selectional transparency in magnetic field. 

This is a novel effect which, as we hope, can be observed at comparatively 
weak magnetic field B < Bj = (po/'jd 2 parallel-to- layers. It was shown 
experimentally ||361 , [27fl, fl28j that the pinning of parallel vortices is much 
weaker than the pinning of pancakes. Therefore, it can be expected that the 
parallel vortices form a highly regular lattice even at fields B < Bj. As it 
was discussed earlier, the PF is not influenced much by the magnetic field in 
this region. An incident EM wave is modulated by the vortex lattice. The 
modulated wave has the same frequency as the incident one, but its wave- 
vector coincides with one of the vortex lattice reciprocal vectors. A most 
important situation is that when the reciprocal vector is parallel to layers. 
Neglecting the in-plane anisotropy, an in-plane wave- vector is determined by 
an integer k and has the length q = 2nk/a. Thus the resonance is shifted to 
a frequency; 

u = ^l l + avW (30) 

where a is a constant of the order of unity. The resonance frequency can be 
regulated by a variation of the magnetic field. If the frequency is fixed, the 
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centers of transparency windows correspond to following values of fields: 



Sfc " {2nv F akf (31) 

The transparency can be observed if the electric vector has a component 
perpendicular to the layers. This phenomenon can be used for the diagnostics 
of the vortex lattice. For example, if the vortex lattice is pinned by the 
periodic potential of the crystal lattice |37j, [j38fl , then, instead of Eq. (|3l"l) , 



the resonance frequency obeys a following equation: 



where So is a value of the field at which the vortex lattice has been quenched. 
The effect can be used for creating magnetic field regulated infrared and 
micro-vawe optical gates. 



10 Conclusions. 

1. In layered structures the PF is strongly reduced, if electric field is perpen- 
dicular to the layers. 

2. The plasma oscillations in SC can propagate even if they are overdamped 
in the normal state. 

3. The experiments cited in the text convincingly show the existence of the 
plasma edge and the resonance plasma attenuation in three HTSC. 

4. The plasma oscillations can propagate in rather clean SC with nodes of 
the OP, since the Cooper-pair breaking is dynamically forbidden. 

5. In Josephson-linked SC, sufficiently high magnetic fields reduce the PF 
dramatically. 

6. At weaker fields, parallel to the layers, an effect of the selectional trans- 
parency should be observable. 

7. Measurement of the magnetic field dependence of the PF can be used for 
the diagnostics of the vortex state in HTSC. 
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Figure Captions 

Fig. 1. 

Theoretical curves for the reflectivity vs frequency measured in units Q p i |J. 
The high-frequency dielectric constant eoo = 25. Different curves correspond 
to different values of the scattering rate: Q p it = 0.01; 0.1; 0.5; ; 1.0; 5.0. 
Fig. 2. 

The plasma edge in YBCO @. Far infrared reflectivity at 300 K (solid 
line), 110 K (dashed line), 60 K (dotted line), and 10 K (dash-dotted line). 
Upper panel: electric field parallel to the a6-plane. Lower panel: electric field 
parallel to the c-axis. 
Fig. 3. 

The plasma edge in La 2 - x Sr x CuOi [[J. The electric field is parallel to the 
c-axis. The values of the concentration x and temperature T are indicated 
on the graphs. 
Fig. 4. 

The resonance attenuation in _Bi 2 S , r 2 CaCM208+<5 H- The surface resistance 
R s vs magnetic field observed at 30 GHz and 2.8 K. A remarkable hysteresis 
is caused by flux trapping in the crystal. In the inset the geometry of the 
experiment is shown. 
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